Twist instability in strongly correlated carbon nanotubes 
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We show that strong Luttinger correlations of the electron liquid in armchair carbon nanotubes 
lead to a significant enhancement of the onset temperature of the putative twist Peierls instability. 
The instability results in a spontaneous uniform twist deformation of the lattice at low temperatures, 
and a gapped ground state. Depending on values of the coupling constants the umklapp electron 
scattering processes can assist or compete with the twist instability. In case of the competition the 
umklapp processes win in wide tubes. In narrow tubes the outcome of the competition depends 
on the relative strength of the e-e and e-ph backscattering. Our estimates show that the twist 
instability may be realized in free standing (5, 5) tubes. 
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The band structure of armchair carbon nanotubes is 
metallic [1]. However, due to the Peierls mechanism 
the metallic state may become unstable at low temper- 
atures. For undoped armchair tubes the possible in- 
stability modes are a lattice deformation with a finite 
wave vector [2] and the twist instability at zero wave 
vector [3l H]. Previous studies of this instability were 
carried out either in the noninteracting electron approxi- 
mation [51 131 m [5J |B] , yielding a very low instability tem- 
perature < OAK even for the thin (5,5) nanotubes, or 
using ab initio calculations in the local density approx- 
imation [7]. These treatments can not account for the 
Luttinger liquid effects, which are expected to be pro- 
nounced in these systems [SJ [^ . 

Here we study the possibility of the Peierls instability 
in armchair carbon nanotubes in the presence of electron- 
electron (e-e) interactions. We find that in the absence of 
backscattering processes the Luttinger liquid correlation 
of electrons significantly enhance the onset temperature 
for the twist Peierls instability. The enhanced transi- 
tion temperature in narrow tubes can be easily reached 
experimentally. We show that depending on the values 
of the e-e umklapp coupling constants the umklapp pro- 
cesses [5J [HI Eni E] either enhance or compete with the 
twist instability. The energy gap due to the umklapp 
processes increases slower with decreasing tube radius R 
than the gap due to the Peierls instability. In case of 
the competition the twist instability can occur only in 
narrow tubes. Based on the theoretical values [T5J [T3] 
of the electron-phonon couplings we estimate the Peierls 
gap in (5,5) tubes with R ~ 0.35 nm to be of the order 
of ~ 70 K. The gaps due to the umklapp processes are 
much harder to evaluate theoretically because of their 



strong sensitivity to the short range part of the Coulomb 
repulsion potential jTOl [TT]. Gate voltage scans of con- 
ductance for armchair tubes [TS] [TH] indicate that the 
umklapp gaps for (5, 5) tubes are of order or below our 
estimates for the twist gap. This suggests that the Peierls 
instability may be realized in torsional nanomechanical 
resonators |17j based on (5,5) armchair tubes. 

We consider an undoped free standing {N, N) arm- 
chair carbon nanotube, see Fig. fl] (a). It is symmetric 
with respect to reflection in the plane that goes through 
the tube axis and maps the A and B sublattices onto 
each other. Thus both electron and phonon modes are 
characterized by parity, ±1 [3]. The low energy electron 
spectrum is formed by two bands with opposite parity 
and zero angular momentum along the tube axis. These 
bands intersect at the two Dirac points as shown in Fig.fl] 
(b). 

Only the phonon modes with zero angular momentum 
along the tube axis couple to the low energy electron 
modes. Furthermore, only the backscattering part of the 
electron-phonon (e-ph) coupling can lead to the Peierls 
instability and is strongly enhanced by the Luttinger liq- 
uid correlations. We therefore neglect the forward scat- 
tering part of the e-ph coupling [5^. The backscatter- 
ing by positive parity phonons involves large momen- 
tum transfer and corresponds to phonons with rather 
high frequencies, ~ lOOOii'. This cuts off Luttinger liq- 
uid renormalization and results only in a finite soften- 
ing of these phonons [TS]. The Peierls instability arises 
from the coupling to the negative parity phonons with 
small momentum. There are three such phonon modes: 
transverse acoustic (TA) and longitudunal/radial opti- 
cal (LO/RO) [3J|T3j. The e-ph coupling arises from the 




FIG. 1: (a) Schematic picture of an armchair tube. FiUed 
and open circles denote atoms in the A and B sublattices. 
The arrows show atomic displacements in the TA and LO 
phonon modes, (b) Free electron spectrum of armchair carbon 
nanotubes: + and — denote the parity of the bands, a — ±1 
denotes the two valleys with the Fermi points at ikp. 



change in the hopping matrix element due to the atomic 
displacements. The coupling to the RO mode is due to 
the curvature of the tube and is smaller than the cou- 
pling to the LO and TA modes in 1/7V [13]. We neglect 
it below. Its inclusion would only enhance the twist in- 
stability. 

The phonon Hamiltonian can be written as 
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^ (K(g)|2 + c.,(gfK(g)|2). (1) 



=L.T 



Here p is the mass per unit length of the tube and Ua{q) 
are Fourier components of atomic displacements. The in- 
dex a — L corresponds to the LO mode with momentum 
independent frequency ujLiq) — uJq and a = T corre- 
sponds to the TA mode with frequency aJriq) — stQ, st 
being the speed of sound. 

The free electron Hamiltonian can be written as 



Hoe 



-ihvp 2^ / 



dx%p'^^„{x)rdxi}ara{x). (2) 



Here Wi? sa 8 x 10^ ra/ s is the Fermi velocity, a — ±1 is 
the valley index, r — ±1 represents right and left movers, 
and a is the electron spin. 

The forward scattering part of the e-e interaction is 
much stronger than the backscattering part, the latter 
being small in 1/A^ [19]. The former can be written as 

ara 

Here we assumed that the Coulomb interaction is 
screened by e.g. a gate at a distance d, and the coupling 
constant being gp « 2e^ ln{d/R). 

The above electron Hamiltonian is bosonized by the 
standard procedure (see, for example, [20jl2T]) 



Yara y^) 



F„ 



exp{-iV^[ea^(a;) - r <^a„{x)]}. (4) 



Here F^a are the Klein factors, ^ is a short distance of 
order of the radius of the nanotube, and $ and 6 are 
boson fields satisfying the commutation relation 

[<^aa{x),Qa'a'{x')] = -iSaa' Saa'0{x - x') , (5) 

with 9(x—x') being the step function. Introducing charge 
and spin modes combining different valleys, 

^aa = [*c+ + a^c- + o-*s+ + aa^s-] /2, 
0a<T = [ec+ + aOc-+(^Qs+ + a<yQs-]/2, 

we rewrite the forward scattering part of the electron 
Hamiltonian as 

H„,+Hp = Y,^ I dx[K-\v^,r+K,ive,n 

3 

(6) 
Here j ~ c±,s±, Uj ~ vp/Kj with Kj being the Lut- 
tinger parameter is the velocity of the j'-th mode. For 
the three modes j = c—,s±, Kj = 1. For the charge 
mode ifc+ = I/a/I + Agp/irvph. Below we will use the 
high charge stiffness approximation, Kc+ ^ 1. 

The coupling of electrons to the TA and LO phonons 
is described by the Hamiltonian 



H, 



dxM{x)[gT'^UT{x) + gLUL{x)], 



(7) 



with M{x) = -iEar^"^V'L^(a;)V'a- 
form of the operator M is 



The bosonized 



M 
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Jl cos(0r$cz.)sin(V7r$si/) 



y=± 



Jl sin(^/7f$cl.)cos(^/7f$s^) 



v=± 



(8) 



By writing M in this form we adopted the convention 
Fi„Faa = 1. For Kc+ < 1 the c-ph coupling ^ 
leads to a strong renormalization of the phonon mode. 
To the second order in the c-ph coupling the phonon 
propagator matrix becomes I?(w„,g) = \Dq {ujn,q) + 
T,{ujn,Q)]~^ , where ujn — 2mrT is the Matsubara fre- 
quency, VQ'^{ujn,q) = da,a'P['^l +ujl{q)] is the bare 
phonon propagator, and T,{ujn,q) is the self energy due 
to the e-ph interaction. The latter has the form 



S(t^„,g) = -P{uJn,q) 
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WLQrq 



-WT9Lq 



(9) 



where P(w„,g) = /dre*'^"^"*?^(M(0)M(r)) and (...) 
denotes thermal averaging with respect to the forward 
scattering Hamiltonian ([6| . In the bosonized representa- 
tion it is given by 



P{^n,q) 



2d^r 



^zuj„T-iqx- f E, ([*j ir)--ij (0)]") 



In the long wavelength limit, q,uj -^ and for Kc+ ^ 1 
we obtain, 



^ fPypy 1 

Vp \ Tr£^ J TT^ 



P{0,0) = P = -^ [ '^ ] ^52(3/8,1/4), (10) 



where B{a, b) is the Euler Beta function. 

The renormalized phonon frequencies are given by the 
poles of det (I>(w„, q)) analytically continued to real fre- 
quencies. The instability first appears when the renor- 
malized frequency of the acoustic mode, ljt turns to zero. 



,~,2 



1 - glvpP 



0, 



(11) 



where gx = 9tI \J ps\vp and g^ = g^/ ^J pujIvF are the 
characteristic dimensionless e-ph coupling constants for 
TA and LO phonons respectively. The mean field twist 
instability temperature is then 



Tc = (5t + 9l) -^B 



(12) 



In one dimension fluctuations shift mean field instabil- 
ities to zero temperature. In the given case the order 
parameter Vu is real and the instability is of the Ising 
type. As is the case for the Ising model, a finite or- 
der parameter exists only at T = whereas at finite T 
we have a state with a finite density of solitons (domain 
walls separating areas with different sign of Vu). The 
mean field Tc is, however, not devoid of meaning and 
gives an estimate of the spectral gaps associated with the 
strong electron-phonon coupling. The exact composition 
of the excitation spectrum is currently unknown. Using 
the Feynmann's variational principle we can estimate the 
magnitudes of the gaps of the singlet excitations of fields 
$j [20]. These excitations represent small fluctuations 
around the minima of the action. In addition to them 
there are kinks of these fields. More complete analysis 
of the spectrum will be given in the extended version of 
this paper. For the condensed state, we assume the lat- 
tice has a static deformation, Vmt — Vt '^l — Ci and 
for the electronic degrees of freedom we use a variational 
action of the form. 
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dxdr 
KjUj 



$1 



«?(V$)^ 
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(13) 



where Aj are variational parameters. The variational 
free energy for electrons is F — F^ + T{S — So)o, where 
(. . .)o denotes averaging with respect to the trial action 
( 13 1. Minimizing the total free energy of both the lattice 



and electrons with respect to Aj and the static lattice 
deformations C and rj one obtains the following equations 



in the Kc+ — > limit. 



1 



_AgT T-r / Aj 






(14) 
(15) 



A? 



1 



2i^ - ^(.9T^/+5^o n 



vpA 



(16) 



This gives identical gaps for the four electron modes. 



^.^YL ( ^(9T + 9l) 

and the spontaneous twist angle 

gr A 



1] 



^(52+51)7^4^ 



(17) 



(18) 



In the framework of the tight binding model, the values 
of the coupling constants can be expressed [121 ESI 1221 [22] 
in terms of the derivative, qJ^ , of the transfer integral 
J{r) with respect to the bond length r. For armchair 
nanotubes, these coupling constants were obtained to be: 



9t 



V3 



„ 9J{r) 
4 " Or 



and gL 



3 9J{r) 
2 dr ' 



where a = 2.5 A is 



the graphene lattice constant and q^ — — AJq/c with 
Jo — 2.6eV the hopping integral, and c = 1.4 A the bond 
length [23]. Here A is a dimensionless constant, whose 
theoretical value is 2 [TJ]. The linear mass density for 
an {N, N) armchair nanotube is p = ANM/a with M 
being the carbon atom mass. The twiston phonon ve- 
locity St ^ 1.4 X lO^'m/s and the LO phonon energy 
ijjQ ^ 0.18 eF [3]. Thus both gx and cjl are propor- 
tional to 1/y/R. For a (5,5) armchair nanotube with 
R ^ 0.35 nm, the transition temperature from Eq. (12 1 



is about 40 i^ and the gap from Eq. (17 1 is about 70 K 



As expected, the mean field instability temperature and 
the gap have the same order of magnitude. The twist 
angle is 77 ~ 3 x 10~^. 

As mentioned above the soliton effects will restore the 
symmetry at finite temperatures. However, if the length 
of the tube is shorter than the typical intersoliton dis- 
tance, at experimental time scales the system will appear 
twisted in one direction. The estimated characteristic 
twist angle for (5,5) nanotubes, ~ 0.01°, is too small to 
be detected by STM imaging. However, for a freely sus- 
pended tube clamped at one end the accumulated rota- 
tion angle at a distance x from the clamp, 9{x) — r]x/R, 
becomes substantial for x ~ 1 fim. This twist can be 
detected in torsional nanomechanical resonators similar 
to the those studied in Ref. [Ill by measuring a defiection 
angle of a paddle attached to an armchair tube. 

The discussion above ignored the backscattering part 
of the e-e interactions. The most relevant backscatter- 
ing terms correspond to the so-called umklapp processes. 



which transfer two right-moving electrons into left mov- 
ing ones or vice versa [3 |9j [lOj [11] . The bosonized form 
of the umklapp interaction can be written as [11] 

Hu = - ^, ^.2 / '^2^cos(\/47r$c+) 

|g3Cos(V4^e,_) + (.93 - 5i) cos(V4^$s+) 
-f5i[cos(V4^$c-) -cos(V4^$s-)]} , (19) 

where the coupling constants can be expressed in terms 
of the matrix elements of the Coulomb interaction of elec- 
trons on the two sublattices with a large {q ^ 2A:f) and 
small {q ~ 0) momentum transfer: g^ = VAA{'^kp) + 
VAB{'2kp) and 31 = Vaa{Q) - Vab{Q) [H [II] . In 
the absence of electron phonon coupling the umklapp 
processes lead to an insulating state with four differ- 
ent gaps for the different modes [H]. The properties 
of the ground state and the magnitudes of the gaps de- 
pend on the values of the coupling constants gi and 173. 
These constants are very difficult to evaluate because 
they depend very sensitively on the short range part of 
the Coulomb interaction and on the spatial dependence 
of the electron density in the atomic orbitals. Within 
the model of point-like atomic orbitals and Coulomb in- 
teraction with a short distance cutoff numerical evalu- 
ation [I^ of the matrix elements gives 33 > 51 > 0. 
In this case it is easy to see from Eqs. ([T9| and ((sl) 



that the umklapp processes compete with the twist in- 
stability. Indeed the umklapp interaction favors the 
condensation of the fields (<i>c+, $c-, *I's+j ©s-) in the 
ground state to (0,0,0,0), while the twist instability re- 
sults in condensation of fields (<i>c+, <&c-,*&s+, *&s-) to 
(0,0,^,^) or (^,^,0,0). We note that for 51 > 
53 > the umklapp processes enhance the twist insta- 
bility and favor the ground state condensation pattern 
(<I>,+ ,<i>,_,<i>,+ ,<l>._) = (0,0,^,^). 

For 53 > gi > the fate of the low temperature state 
thus depends on the relative strength of the e-ph and 
umklapp coupling constants. To estimate which ground 
state has the lowest energy one can compare the umklapp 
gap with that due to the twist instability. For A'c+ ^ 1 

the former scale as A„ ^ ^^ whereas the latter as A ~ 

2 
^^^Sr- Therefore in wide tubes the umklapp processes 

win. The outcome of the competition in narrow tubes 

depends on the values of the umklapp coupling constants. 

Theoretical estimates [10] for the gaps range from tens 

to hundreds of Kelvins for (5,5) tubes. Experimentally, 

armchair carbon nanotubes with R ~ 0.8 nm show no 

sign of the gap down to T = AK [TB]. Using l/B? scaHng 

we infer that the umklapp gaps in (5,5) armchair tubes 

should be below our estimates of 70 K for the twist gap. 

These estimates suggest that twist instability may occur 

in freely suspended (5, 5) tubes. It could be observed 

by measuring the twist angle in the torsion resonators 

similar to those studied in Ref . HT] 
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